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Abstract. We investigate the effect of the anisotropy of a harmonic trap on 
the behaviour of a fast rotating Bose-Einstein condensate. This is done in the 
framework of the 2D Gross-Pitaevskii equation and requires a symplectic reduction 
of the quadratic form defining the energy. This reduction allows us to simplify 
the energy on a Bargmann space and study the asymptotics of large rotational 
velocity. We characterize two regimes of velocity and anisotropy; in the first one 
where the behaviour is similar to the isotropic case, we construct an upper bound: 
a hexagonal Abrikosov lattice of vortices, with an inverted parabola profile. The 
second regime deals with very large velocities, a case in which we prove that the 
ground state does not display vortices in the bulk, with a ID limiting problem. In 
that case, we show that the coarse grained atomic density behaves like an inverted 
parabola with large radius in the deconfined direction but keeps a fixed profile 
given by a Gaussian in the other direction. The features of this second regime 
appear as new phenomena. 
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1. Introduction 

Bose-Einstein condensates (BEC) are a new phase of matter where various as- 
pects of macroscopic quantum physics can be studied. Many experimental and 
theoretical works have emerged in the past ten years. We refer to the monographs 
by C. J.Pethick-H.Smith [TTj, L.Pitaevskii-S.Stringari [18] for more details on the 
physics and to A.Aftalion [2] for the mathematical aspects. Our work is motivated 
by experiments in the group of J.Dalibard [11] on rotating condensates: when a 
condensate is rotated at a sufficiently large velocity, a superfluid behaviour is de- 
tected with the observation of quantized vortices. These vortices arrange themselves 
on a lattice, similar to Abrikosov lattices in superconductors [T]. This fast rotation 
regime is of interest for its analogy with Quantum Hall physics [SI El EI] • 

In a previous work, A.Aftalion, X.Blanc and F.Nier [3] have addressed the math- 
ematical aspects of fast rotating condensates in harmonic isotropic traps and gave 
a mathematical description of the observed vortex lattice. This was done through 
the minimization of the Gross-Pitaevskii energy and the introduction of Bargmann 
spaces to describe the lowest Landau level sets of states. Nevertheless, the ex- 
perimental device leading to the realization of a rotating condensate requires an 
anisotropy of the trap holding the atoms, which was not taken into account in [3] . 
Several physics papers have addressed the behaviour of anisotropic condensates un- 
der rotation and its similarity or differences with isotropic traps. We refer the reader 
to the paper by A. Fetter [8], and to the related works [Ml [19], [20]. The aim of the 
present article is to analyze the effect of anisotropy on the energy minimization and 
the vortex pattern, and in particular to derive a mathematical study of some of 
Fetter's computations and conjectures. Two different situations emerge according 
to the values of the parameters: in one case, the behaviour is similar to the isotropic 
case with a triangular vortex lattice; in the other case, for very large velocities, 
we have found a new regime where there are no vortices, and a full mathematical 
analysis can be performed, reducing the minimization to a ID problem. The exis- 
tence of this new regime was apparently not predicted in the physics literature. This 
feature relies on the analysis of the bottom of the spectrum of a specific operator 
whose positive lower bound prevents the condensate from shrinking in one direction, 
contradicting some heuristic explanations present in [8]. Our analysis is based on 
the symplectic reduction of the quadratic form defining the Hamiltonian (inspired 
by the computations of Fetter [8]), the characterization of a lowest Landau level 
adapted to the anisotropy and finally the study of the reduced energy in this space. 

1.1. The physics problem and its mathematical formulation. Our problem 
comes from the study of the 3D Gross-Pitaevskii energy functional for a fast rotating 
Bose-Einstein condensate with particles of mass m given by 

(1-1) ^GP(0) = W, 0)l2(m3) + ^^UWli^M^), 

where the operator Ti. is 

{1.2) n = ^{h'Dl + h'Dl + h'Dl) + ^ {u^lx' + + uy) - QixhD, - yhD^), 

where h is the Planck constant, = {2iTT)~^dx, ujj is the frequency along the j-axis, 
fl is the rotational velocity, and the couphng constant g^d is a positive parameter. 
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In the particular case where oo^ = ojy^ the fast rotation regime corresponds to 
the case where VL tends to and the condensate expands in the transverse di- 
rection. It has been proved [1] that the minimizer can be described at leading 
order by a 2D function tlj{x,y), multiplied by the ground state of the harmonic 
oscillator in the 2;-direction (the operator h"^ / {2m)Dl + mu;^z^/2), which is equal 
to {2muJzh~^Y/^e~'^^'^^'^ ^ . This property is still true in the anisotropic case if 
ujy <^ uj^. The reduced 2D energy to study is thus 

(1-3) ^{'IP) = (W,V^)l2(R2) + ^y^||^||i4(]R2), 

where the operator Ho is 

1 777 

(1.4) Ho = ^(h'Dl + h'Dl) + -[ujIx^ + uoy) - n{xhDy - yhD,), 

and the coupling constant g2d takes into account the integral of the ground state in 
the 2;-direction: 

(1.5) g2dN = , where g is dimensionless (and > 0). 

772 

Since h has the dimension energy x time, it is consistent to assume that the wave 
function ip has the dimension 1/length, with the normalization ||'?/'||i2(]K2) = 1. We 
define the mean square oscillator frequency uj± by 

2 ^ / 2 , 2\ 

^± = ^y'^x + ^y) 

and the function u by 

(1.6) ^Ij{x, y) = h~^'^m^'^uj'l\{h-^'^m"^uj'l\, h-^'^m^'^u'l^y), 
so that 

lkllL2(iR2) = ||^/'||l2(r2) = 1, 92dN\\i)\\'\^A(^2) = ghu±\\u\\l4^^2y 

We also note that the dimension of h^^^'^m^^'^uj]^'^ is 1/length, so that 

Xi = h~^^'^m^^'^ujY'^x, X2 = h^^^'^m^^'^ujY'^y^ u{xi,X2) are dimensionless. 
Assuming u"^ < Uy, we use the dimensionless parameter u to write 

2/1 2\ 2 2 /I I 2\ 2 

and we get immediately 

■'■ lii" i|2 ,-'-|ln„.ll2 .-'-/I ,,2mi„.i|2 |-'-/i,,,2mi„„.i|2 



= ^PiM||i2(R2)+-p2M|L2(K2)+-(l-z/ )||a;iM|L2(iR2)+-(l+i/ )||a;2M||i2(R2) 



hu^ 2" ^ "^^^^ 2 



Finally, we have 



- —{iXiD2 - X2-Di)m,m)l2(r2) + ^||M||i4(]R2). 



(1.7) 7— •= Egp{u) = {Hu, u) + hu\\i, 

(1.8) 2H = Dl + Dl + {l-u^)xl + {l + u'^)xl-2u{xiD2-X2Di), u = —. 

UJ± 
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where uj,i',u,g are all dimensionless and 2:2(-ir2) = 1. The minimization of this 
functional is the mathematical problem that we address in this paper. The Euler- 
Lagrange equation for the minimization of Egp{u), under the constraint ||m||x,2(]k2) = 
1, is 

(1.9) Hu + g\u\'^u = Xu, 

where A is the Lagrange multiplier. We shall always assume that < u"^, i.e. 
uj"^ + < 1 and define the dimensionless parameter e by 

(1.10) u;^ + u^ + e^ = l. 

The fast rotation regime occurs when the ratio ^"^/uj^ tends to 1_, i.e. e tends to 0. 
Summarizing and reformulating our reduction, we have 

I g r 

(1-11) Egp{u) = -{q^,u,e^,u)L2(M.2) + - / |m|^c?x, 

where quj,u,e is the quadratic form 

(1.12) g^,,,,(xi, X2, ^1, 6) = ^1 + ^2 + (1 - ^^Wi + (1 + - 2c^(2;i6 - X2^i), 

which depends on the real parameters u^u^e such tha^ (11.101) holds. Here q^^^ is 
the operator with Weyl symbol qu],v,ei that is: 

(1.13) = Dl + Dl + (1 - v^)xl + (1 + v'')xl - 2uj{xiD2 - X2D1), 

where Dj = dj/{2i7r). We would like to minimize the energy Egp{u) under the 
constraint \\u\\l2 = 1 and understand what is happening when e ^ 0. 

1.2. The isotropic Lowest Landau Level. When the harmonic trap is isotropic, 
i.e. when = 0, it turns out that, since u"^ + e"^ = 1, 

(1.14) q = g^,o,e = (6 + ^2:2)' + (6 - u}Xif + e\xl + xl) 
so that 

1 uo e"^ 

Egp{u) = -\\{Di + ujX2)iJ + i{D2 - ujxi)uf + —\\uf + — 

I In A ^ 

We note that, with z = xi + 1x2-, 

Dl + UJX2 + i(D2 — uJXi) = — d — iujz = — (B + nuz), 

ZTT 111 

hence the first term of the energy is minimized (and equal to 0) if m G LLL^-i, 
where 







' 1 !■ 


a; n p H 




\u\ 









;i.l5) LLL^-i = {ue L'\M:'),u{x) = /(^)e-""'"''} = ker(a + ttujz) n L^^^^ 



with / holomorphic. We expect the condensate to have a large expansion, hence the 
term J to be small. Thus, it is natural to minimize the energy Egp in LLL^-i. 
It has been proved in [1] that the restriction to LLL^^i is a good approximation 



""^Of course there is no loss of generality assuming that e, ly are nonnegative parameters; we may 
also assume that oj > 0, since the change of function u{xi, X2) ^ u{~xi, X2) preserves the L^-norm, 
is unitary in L^, corresponding to the symplectic transformation {xi, X2, Cii C2) ^ X2, — Cii C2) 

and leads to the same problem where uj is replaced by —uj. 
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of the original problem, i.e. the minimization of Ecp in L^(M^). We get for u G 
LLL^-i, \\u\\l2 = 1, 

1 uj e'^ 

Egp{u) = -\\ {Di + ujX2)u + i{D2 - ujxi)u f + — + — | 
z ^ Ztx I 

{in) ^ ^ {d+nijjz)u=0 

and with u{x) = v[{ujeY^'^x){ujeY^'^ (unitary change in L^(]R^)), 

^ ^ ' ' ' |2| / m2 T , 2 / I / m4 



^gp(m) = ^ + ^ |l/|>(l/)|% + ^'^7 J \v{y)rdyj. 

The minimization problem of Ecpiu) in the space LLL^-i is thus reduced to study 

(1.16) Elll{v) = II \x\v\\l, + uj^g\\v\\l,, v e LLL„ 

i.e. with 2; = xi + ix2, v{xi,X2) = f{z)e~'^'^ / entire (and v G L^(M^)). This 
program has been carried out in the paper [3] by A. Aftalion, X. Blanc, F. Nier. 
In the isotropic key point is the fact that the symplectic diagonalisation of 

the quadratic Hamiltonian is rather simple: in fact revisiting the formula (I1.14p . we 
obtain easily 



+ (^)(ei + x2r + (^)(6-xi 




= 2-1/2(1 _^)l/2(^^_ ^2), m = l-'^, 2/1= 2-1/2(1 _^)-l/2(^2+^^)^ 
= 2-1/2(1+^)1/2(^^+^2), fi2 = l+u:, 2/2 = 2-1/2(1 + ^)-V2(^^_ ^2), 

SO that the linear forms {yi,y2^ViiV2) are symplectic coordinates in R"^, i.e. 

{^1, z/i} = {V2, z/2} = 1, {^1, ^2} = {viiy2} = {v2, yi} = {yi,y2} = o. 

In [3], an upper bound for the energy is constructed with a test function which 
is also an "almost" solution to the Euler-Lagrange equation corresponding to the 
minimization of fll.l6|] in LLL^. This almost solution displays a triangular vortex 
lattice in a central region of the condensate and is constructed using a Jacobi Theta 
function, which is modulated by an inverted parabola profile and projected onto 
ILL,. 



1.3. Sketch of some preliminary reductions in the anisotropic case. The 

analysis of the reduced energy in the anisotropic case yields two different situations: 
one is similar to the isotropic case and the other one is quite different, without 
vortices. To tackle the non-isotropic case where > in (11.131) . one would like to 
determine a space playing the role of the LLL and taking into account the anisotropy. 
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Step 1. Symplectic reduction of the quadratic form q^^y^^. Given the quadratic form 
(lu),v,e (ll-l^P - identified with a 4 x 4 symmetric matrix, we define its fundamental 
matrix by the identity F = —cr~^q^^y^^ = aq^^^^e where 

cr = ^ ^ "1^^^ is the symplectic matrix given in 2 x 2 blocks. 

The properties of the eigenvalues and eigenvectors of F allow to find a symplectic 
reduction for ga;,i^,e- 

Step 2. Determination of the anisotropic LLL. The anisotropic equivalent of the 
LLL can be determined explicitely, thanks to the results of the first step. We find 
that it is the subspace of functions u of L^(]R^) such that 



(2 2 \ 

- ^) + {^2X2)^(1 + ^)] j exp 



-I— — X1X2) 
4a 



where / is entire. The positive parameters a,7,/32 are defined in the text and 
are explicitely known in terms of a;, z/. We also determine an operator M, which 
can be used to give an explicit expression for the isomorphism between L^(R) and 
the anisotropic LLL as well as to express the Gross-Pitaevskii energy in the new 
symplectic coordinates. 

Step 3. Reseating. Introducing a new set of parameters (cj, z/, e are positive satisfying 
(HlOj), g>0 given by (O)), 

{1.19) nl = {2u' + e'){l + ^^), a = V^AT4i7^ gi = 9^^^f^, 

1.20 K = go = = —,(32 = — , ^ 2 I 2 ' /i2 = l + w+tt, 
jj2 4p2 uj a + ziv^ + 

we show that, after some rescaling, the minimization of the full energy Egp{u) of 
(11.111) can be reduced to the minimization of 

(1.21) E{u)= f l{e'xi + K-'xl)\u\'' + ^\u\\ 



on the space 

(1.22) Aq = {ue L2(M^), u{xi,X2) = /(2)e-^l^l'/^ / holomorphic, z = Xi+ix2}. 

The point is that, after some scaling, we are able to come back to an isotropic space. 
The orthogonal projection IIo of L^(]R^) onto Aq is explicit and simple: 

(1.23) (nou)(x) = [ e^^l"^^l'+'"("^^^-^^"^)u(?/)di/. 



We are thus reduced to the following problem: with E{u) given by (11.211) . study 
(1.24) J(e, k) = inf {E(m), u E Aq, \\u\\L^m2) = l}. 

The minimization of E without the holomorphy constraint yields 

1.25 M = — — 1 - -2 - -2 , where i?i = — -] , R2 = [ — j] ■ 
nRiR2 Rf R2 V / V '^^ J 

As e tends to 0, i?i always tends to infinity (in fact Ri > e^^^^), but the behaviour 
of i?2 depends on the respective values of e and k, that is of e and u. 
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Step 4- Sorting out the various regimes. Recalling that the positive parameter ly 
stands for the anisotropy, we find two regimes: 

• z/ <^ e"*^/^ (weak anisotropy): R2 ^ oo (in fact, i?2^^ ~ min{e~'^^^,e^^^i'^^)). 
Numerical simulations (Figure [T]) show a triangular vortex lattice. The behaviour 
is similar to the isotropic case except that the inverted parabola profile f ll.25p takes 
into account the anisotropy. We will construct an approximate minimizer. 
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Figure 1 . Plot of the zeroes of the minimizer (left) and the density (right) 
for = 0.002, v = 0.03. Triangular vortex lattice in an anisotropic trap. 



• u ^ (strong anisotropy): R2 (in fact R2 ~ e^^^u"^). Numerical 
simulations (Figure ^ show that there are no vortices in the bulk, the behaviour 
is an inverted parabola in the Xi direction and a fixed Gaussian in the X2 direction. 
Thus, the size of the condensate does not shrink in the X2 direction and fll.25p is not 
a good approximation of the minimizer. The shrinking of the condensate in the X2 
direction is not allowed in Aq (see (11.221) ) because the operator X2 is bounded from 
below in that space by a positive constant and the first eigenfunction is a Gaussian 
in the X2 direction. We find an asymptotic ID problem (upper and lower bounds 
match) which yields a separation of variables. 
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Figure 2. Plot of the zeroes of the minimizer (left) and the density (right) 
for = 0.002, v = 0.73. No vortex in the visible region. 



1.4. Main results. 

1.4.1. Weakly anisotropic case. In a first stepH, we assume tliat, witli k given by 
(IL20D, 

(1.26) e < K < e^'^. 

Tlie isotropic case is recovered by assuming k = e. This case is similar to the isotropic 
case and we derive similar results to the paper [3], namely an upper bound given by 
the Theta function but we lack a good lower bound. 

We recall that the Jacobi Theta function 0(z, r) associated to a lattice Z © Zr is 
a holomorphic function which vanishes exactly once in any lattice cell and is defined 
by 

(1.27) e(z,r) = - ^ (-l)V("+^/2)2e(2"+i)-'^, zeC. 

n=— oo 

This function allows us to construct a periodic function on the same lattice: Ur is 
defined by 

(1.28) M^(xi,X2) = e^('''"l'^l')e (v/T72;,r) , z = xi+ix2, r = rR + iTi, 
\Ur\ is periodic over the lattice Z © rZ, and Ur satisfies 

(1.29) Ho = A^u^, 
with 



;i.30) 



A. 



^We shall see that k v + s in the sense that the ratio k/{v + e) is bounded above and below 
by some fixed positive constants, so that the weakly anisotropic case is indeed v <C e^^"^. 
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and 

(1.31) 7(r):- ^ 



UK?) 



2 • 



The minimization of 7(r) on all possible r corresponds to the Abrikosov problem. 
It turns out that the properties of the Theta function allow to derive that 

(i,fc)GZ2 

and prove (see [3]) that r i-^ 7(r) is minimized for t = j = e^*'^/'^, which corresponds 
to the hexagonal lattice. The minimum is 

(1.32) 6 = 7(j) ^ 1.1596. 

The function Ur allows us to construct the vortex lattice and we multiply it by the 
proper inverted parabola to get a good upper bound: 

Theorem 1.1. We have for I{e, k) defined in (11.241) . b given in fll.321) . k in fll.201) . 

(1.33) 3 < ^(^' ^) ^ 3 V ^ + 1,^ U 

when (e, ne^^^'^) (0, 0). Moreover, the following function provides the upper bound: 

(1.34) V = Uo{UrP), 
where Ur is defined by f l 1 . 281} with r = and 

'^''^ = ^wi^2\ vm~ J ^'^v^j '""^=1^; ■ 



We expect to be a good approximation of the minimizer and the energy asymp- 
totics to match the right-hand side of (11.331) . Thus, the lower bound is not optimal 
( it does not include b). In addition, the test function (ll.34p (with a general r ^ j 
a priori) gives the upper bound of (11.331) with 7(r) instead of b. The proof is a 
refinement of that in [3]. 

1.4.2. Strong anisotropy. In the case where the rotation is fast enough in the sense 
that 

(1.35) K>£i/3 

we have found a regime unknown by physicists where vortices disappear and the 
problem can be reduced in fact to a ID energy. 

Theorem 1.2. For I{e,K,) defined in (11.241) . b given in (ll.32p . k, in (11.201) . we have 

(1.36) ^ lim (^-^^4/^] = 

where 
(1.37) 

J = inf { j ]fp{tf + I y p{t)\ p real-valued e L^R) n L^(M), Mi^iR) = l}. 
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In addition, if u is a minimizer ofI{e,K), then 

(1.38) ^|n(|^,x,)|->2VV--ip(a;0, 

in L^(]R^) n L^(]R^), where p is the minimizer of J. 
Note that the minimizer p of fll.37p is exphcit: 



A few words about the proof of Theorem ll.2[ The first point is that the operator 
Iloxlno (see (11.221) . fll.23p ) is bounded from below by a positive constant: 

Vn G Ao, / xl\u\^ > — / I^P- 
JR2 47r J^2 

This is proven in Lemma H3] below. Actually, the spectrum of this operator is purely 
continuous, and any Weyl sequence associated with the value 1 / (An) converges (up 
to renormalization) to the function 

(1.39) Uo{xi, X2) = exp {—nxl + z7rxiX2) , 

which satisfies the equation Uq^uq) = -^Uq. This gives the lower bound 

and indicates that in order to be close to this lower bound, a test function should 
be close to fll.391) . Thus, the second point is to construct a test function having the 
same behaviour as fll.39p in X2, and a large extension in Xi. This is done by using 
the function 

«i(xi,X2) = ^e-5^"£e-i((^-^^)'-2^^^^^)p(yi)rfyi, 

which is equal to Ilo{p{xi)6o{x2)) , where 60 is the Dirac delta function and p any 
real-valued function of one variable. This test function is then proved to be close 
to 2^/^e~'^'^2p(xi), which allows to compute its energy, and gives the upper bound, 
provided that p{t) = e^^^p{e'^^^t), where p is the minimizer of f ll.37p . Finally, in 
order to prove the lower bound, we first extract bounds on the minimizer from the 
energy, which allow to pass to the limit in the equation (after rescaling as in (11. 380 ). 
hence prove that the limit is the right-hand side of fll.38p . This uses the fact that 
the energy appearing in fll.37p is strictly convex, hence that any critical point is the 
unique minimizer. 

The paper is organized as follows: in section 2, we review some standard facts on 
positive definite quadratic forms in a symplectic space. This allows us, in section 3, to 
construct a symplectic mapping x, which yields a simplification of the quadratic form 
q. In section 4, quantizing that symplectic mapping in a metaplectic transformation, 
we find the expression of the LLL and manage to reach the reduced form of the 
energy (Proposition 4.5). Section 5 is devoted to the proof of Theorem 1.1 and 
section 6 to Theorem 1.2. 
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Open questions. We have no information on the the intermediate regime where, for 
instance, e^^^/n converges to some constant R^^^ (in that case, Ri ~ e"^/^, R2 ~ Rq). 
We expect that the extension in the X2 direction depends on Rq and wonder whether 
the condensate has a finite number of vortex hnes. We have not determined the 
hmiting problem. 

Acknowledgements. We would like to thank A. L. Fetter and J.Dalibard for very 
useful comments on the physics of the problem. We also acknowledge support from 
the French ministry grant ANR-BLAN-0238, VoLQuan and express our gratitude 
to our colleagues participating to this ANR-project, in particular T.Jolicceur and 
S.Ouvry. 

2. Quadratic Hamiltonians 

We first review some standard facts on positive definite quadratic forms in a 
symplectic space. 

2.1. On positive definite quadratic forms on symplectic spaces. We consider 
the phase space M" x M^, equipped with its canonical symplectic structure: the 
symplectic form a is a bilinear alternate form on M?"^ given by 

(2.1) a((x, 0; (y, r/)) = e ■ y - r/ ■ X = (aX, F), 

(2.) w,th ^"), 

where the form a is identified with the 2n x 2n matrix above given in n x n blocks. 
The symplectic group Spin) (a subgroup of S'/(2n, M)), is defined by the equation 
on the 2n x 2n matrix x, 

(2.3) xVx = cr, i.e. VA, Y G M^", (ax A, xY) = {aX, Y). 

The following lemma is classical (see e.g. the chapter XXI in [10], or [T5]). 

Lemma 2.1. Let B G GL{n.,M.) and let A,C be n x n real symmetric matrices. 
Then the matrix S, given by n x n blocks 

^ _( -B-'C \_fl 0\ (B-^ \ (I -C\ 

[Z.^) B*-AB-^C)~\A B*)\Q I ) 

belongs to Sp{n). Any element of Sp{n) can be written as a product 

N.B. The first statement is easy to verify directly and we shall not use the last statement, 
which is nevertheless an interesting piece of information. For a symplectic mapping H, to 
be of the form above is equivalent to the assumption that the mapping x 1— > pr\r.{x © 0) 
is invertible from M" to W^. 

Given a quadratic form Q on M^"-, identified with a symmetric 2n x 2n matrix, 
we define its fundamental matrix F by the identity 

F = -a~^Q = aQ, so that for A, F G M^" {aY, FX) = {QY, A). 

The following proposition is classical (see e.g. the theorem 21.5.3 in |10j). 
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Proposition 2.2. Let Q be a positive definite quadratic form on the symplectic 
M" X M^ . One can find x ^ Sp{n) such that with 



!>2n 



3X = xY, Y = {yi, . . . ,yn,r]i, . . . ,r]n), 



{QX, X) = {QxY, xY) = Yl + /^^D' > 0- 

The {±i/ij}i<j<„ are the 2n eigenvalues of the fundamental matrix, related to the 
2n eigenvectors {cj ±^£^j}i<j<n- The {Gj,£j}i<j<n rnake a symplectic basis o/M^"; 

a{ej, Ck) = 6j^k, o-{£j, £k) = (^{ej, Ck) = 0, 

and the symplectic planes Uj = Mcj © Wej are orthogonal for Q. 

N.B. A one-line-proof of these classical facts: on C^" equipped with the dot-product given 
by Q, diagonalize the sesquilinear Hermitian form ia. 

2.2. Generating functions. We define on M" x M" the generating function S of 
the symplectic mapping of the form 'E.a^b,c given in the lemma [211] by the identity 

(2.5) S{x,r]) = ^{{Ax,x)+2{Bx,r]) + (Cr/,r/)). 

We have 

In fact, we see directly 

A IJ \ B*J [o I J \ r] J ~ \A IJ \B*r] J ~ \Ax + B*r]J ' 

Given a positive definite quadratic form Q on M."^^, identified with a symmetric 
2n X 2n matrix, we know from the proposition 12.21 that there exists x ^ Sp{n) such 
that 



X*QX = Q , /^^ = diag(/ii, 



Looking for x = '^a,b,c given by a generating function 5* as above, we end-up (using 
the notation q{X) = {QX,X) with X G M^") with the equation 

q{x,dxS) = WfidnSW^ + \\r]\\^, ^dr,S = {fJ'jdn^S)i<j<n e M", 



I" vlCPn 



where || ■ || stands for the standard Euclidean norm on M". This means 
(2.7) q{x,Ax + B*r]) = \\fi{Bx + Cr])f + \\r]f. 

We want now to go back to the study of our quadratic form (11.121) . 
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2.3. Effective diagonalization. 



Lemma 2.3. Let q be the quadratic form on given by fll.l2p . where uj,v,e are 



nonnegative parameters such that uP' 
mental matrix are ±ifii, ±z/i2 with 



e = 1. The eigenvalues of the funda- 



(2.8) 
(2.9) 



< = 1 + Lu'^ - a < fij = 1 + Lu"^ + a, 



a 



2^2 



In the isotropic case v = 0, we recover /ii = 1 — u;, /i2 = 1 + cu. When e > 0, we have 
< < /X2 fl^t^ Q is positive- definite. When e = we have /ii = < ^2, o,nd q is 
positive semi-definite with rank 2 if u = and with rank 3 if u > 0. 



Proof. The matrix Q of g is thus 



(2.10) Q 



\ 



2 

l + 

UJ 









-U) 





CO 

1 





-to 


1 



and F = aQ 



( ° 




1 


0^ 


— OJ 








1 












V 




— w 


0/ 



The characteristic polynomial p of F is easily seen to be even and we calculate 
p{\) = dei{F-\h)=\^ + 2{l+uj^)\^ + {l-uj^f 

The four eigenvalues of F are thus ±2\/ 



± \Jv^ ^ 4a;2, proving the first state- 
ment of the lemma. Since (1 + a;2)2 _ ^ ^ _ _ ^4 ^ e'^{2i''^ + ^2)^ we get 
fif = e2(^2z/2 -|-e2)/^2_ rjj^g statements on the cases u = 0,e > are now obvious. 
When e = = u, we have a; = 1, and rankg = 2 as it is obvious on fll.171) . When 
e = 0, 1/ > 0, we consider the following minor determinant in F, cof actor of /31 

, 

(-l)(-cj^ + + 1) = -2i/V 0, 




,2 



1 

1 

-UJ 



-u- - 1 

so that rankQ = rankF = 3 in that case. 



□ 



N.B. We may note here that the condition uj'^ + u'^ < 1 is an iff condition on the real pa- 
rameters u, UJ for the quadratic form (jl.l2p to be positive semi-definite. This is obvious on 
the expression (jl.lTp in the isotropic case z/ = 0, and more generally, the (non-symplectic) 
decomposition in independent linear forms 

,2\ 



g = (^1 + UJX2f + (6 - wxi)2 + x?(l - 1/2 - u;2) + a;2(l + 1/2 



UJ 



shows that q has exactly one negative eigenvalue when t(;2 + zy2 > 1 > uj"^ — u'^, and exactly 
two negative eigenvalues when u;2 — i/2 > 1. As a result, when ci;2 + ;y2 ^ ^j^g operator 
q^ is unbounded from below. 

Using now the equations (12. 7p . (11.120 and assuming that we may find a linear 
symplectic transformation given by a generating function (12.51) . we have to find 
A, B, C like in the lemma EH] with n = 2, so that for all (x, 77) G IR2 x R2, 

\\Ax + B*r]f +\\x\\'^ + u'^{xl-xl)-2uj{x h{Ax + B*T])) = \\fi{Bx + Cri)f + 

with xAC, = Xi,^2— a^2'Ci! = diag(/ii, ^12). At this point, we see that the previous iden- 
tity forces some relationships between the matrices A, B, C. However, the algebra 
is somewhat complicated and assuming that B is diagonal. A, C are (symmetrical) 
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with zeroes on the diagonal lead to some simplifications and to the following results. 
We introduce first some parameters: 



(2.11) /3i 

(2.12) 132 



a — 2a;2 -)_ y'i 2uj^i 
2ujp.2 a + 2a;^ — v"^ 



a + 2u;2 + y'i 



2uj^2 



2a 

2.13 7 = —, 
(2.14) A? 



since (a - ^w^Y ~y = 4t<J + ^ui - A.u?ol = i^u? ii\ 



since (a + 2up-)'^ - = 4aj^ + 4a;'* + ^u?a = 4tj^/x| 



a - 2a;^ + 



(2.15) Ai 

and we have 
(2.16) 

We define also 



/il + /?l/32M2 1 + 

M2 1 



1 + 



1 



/^2 + /?l/?2/il 1 + 1 + ^ 



/'2 



-2a;g-iy2 
o+2a;2+!y2 



2a 

a + 2u;2 + 

2a ' 



A^ + A2 — 1 H , A^Ag 



(a + 



2\2 



(2.17) d 



7A1A2 



A? + A^ 



2 ' 2A1A2 
Lemma 2.4. M^e define the 2x2 matrices 



which gives cd 



4a2 

2a(l + z/V«) a + 



4a; 



2oj 



5 



A^^ 
Ao- 



rf-i 









Ai A2 



Ai A2 

-cd 



cd 



The 4x4 matrix given with 2x2 blocks by 



X — ^A,B,C 

belongs to Sp{2) and 



h \ fB 



A h 







-C 



(2.18) 



X 



( Ai 







A2 



B* \0 L 







- X2cd CA2 
cAi / 



(2.19) 



/ CA2 





^\ 

^ 

-^ + \icd Ai 



cAi 



\-^ + A2crf A2/ 

Proof. The lemma [2?T] gives that x ^ 'S'p(2) and we have also 

72 C\ /5 



X 



/o / \ 



?* — 1 



I2 

-A /2 



The remaining part of the proof depends on the formula giving 'Ea,b,c in the lemma 
12.11 and a direct computation whose verification is left to the reader. □ 
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Lemma 2.5. Let x be the symplectic matrix given by (12.1 8p and Q be the matrix 
given in fl2.10p . Then, with fij given by fl2.8l) . we have 

(2.20) x*Qx = diag(/i?,/i2,l,l). 



The (tedious) proof of that lemma is given in the appendix 17.3.11 
Using the expression of in (I2.18p . defining 



(yi\ 




2/2 








\V2j 





( 



0X2 








A2\ 

d \ 









cAi 


Ai 
d 







X2 





-i + Aicrf 


Ai 









^ + Aacrf 








X2) 







(2.21) 



we get from the lemma [2751 the following result. 

Lemma 2.6. For (xi, X2, ^1, ^2) £ K^, (z/i, ^2, ''71, ''72) £ g'zi'en 6y (12.2 ip . we /iave 
t/ie following identity, 

fAvl + l^yl +Vl+V2= A^?(cAiX2 + A2d"^C2)^ +/X2(cA2Xi + Aid"^^i)^ 

+ ((-dAg ^ + Aic(i)x2 + AiCi)^ + ((-dA^^ + X2cd)xi + A26)^ 

= + el + (1 - i^')a^? + (1 + - 2uj{xi^2 - X2C1), 

where the parameters c, A2, d, Ai are defined above (note that all these parameters are 
well-defined when [uj, v) are both positive with u'^ + < 1). 



We have achieved an explicit diagonalization of the quadratic form (11.120 and, 
most importantly, that diagonalization is performed via a symplectic mapping. That 
feature will be of particular importance in our next section. Expressing the param- 
eters in terms of a,Lj, p.e (cf. section [77^ . we obtain 

q = (2-^/^a~'l\a - 2u? + v'^fl'^i^ - 2-^/^u;-'a-'l\a - 2lo^ + v^f'\a - v^)x2 
,a^2u? 



(^2-1/2^-1/2 , 



1/=^. 1/2 (21/ V +£4)1/^ 

2z.2 + £2 ) ^2 ^ 



+ 



i^^^!£^±^(aV2 + , 2^-1/2)2-3/2,-1 ^ 

l^2 ^ 2i/2 + e2 



(1 + a;2 + a) V22i/2«-i/2^(a + 2^-" + i/2)-i/2^^ 



+ (1 + ^2 ^ ^)l/2(i ^ ^-1^2)2-l/2„l/2(^ + 2J2 + z.2)-l/2^^ 

+ ( 2-i/2a-V2(a + 2u? + ^.2)1/2^2 - 2-3/2a;-ia-i/2(a - y^){a + 2a;2 + ^2)1/2 



X\ 



so that 
(2.22) 



2a 



2 I o, ,2 



a + 2tu' 



2a;yU2 



^2+ (^T— )^i 



1 2 



+ 2cu^ 



a{a + 2^2 + 1/2) 



a + 2ci;2 + 
2a 



2cj 



■ ^ a - v\ 



'7i 



16 



A. AFTALION, X. BLANC, AND N. LERNER 



The equation (12.221) encapsulates most of our previous work on the diagonahzation 
of q. In the appendix 17.3.21 we provide another way of checking the symplectic 
relationships between the linear forms, UjiTji. 

We have seen in Lemma [2.31 that when e = 0, > 0, the rank of q is 3, whereas 
its symplectic rank is 2. Indeed, e = and > 0, we have 



(2.23) 



2a 

l + uj"^ + a 
a(a + 2cj2 + u^] 



a — V 
2uo 



)X2 



+ 



2 /a + 2tu2 + z/2 



+ 



2a 



2 2 
^2 2/2 



6-( 



a — v 

2uj 



)xi 



3. Quantization 



3.1. The Irving E. Segal formula. Let a be defined on x (say a tempered 
distribution on R^"). Its Weyl quantization is the operator, acting for instance on 



(3.1) 



g2i.(x-x')«a(^±^, ^)u{x')dx'dC 



In fact, the weak formula {a'^u,v) = J^27i(i{x,^)H{u,v){x,^)dxd^ makes sense for 
a E o5^'(R^"), M, f G ^(R") since the Wigner function 7i{u,v) defined by 



n{u,v){x,C) 



e-2*""'«u(a: + |^)tXx - ^)dx' 



belongs to ^(R^") for u,v E ^(R") . Note also our definition of the Fourier 
transform u{C,) = J e~'^^'^^'^u{x)dx (so that u{x) = J e^^'^^'^u{^)d^) and 



1 du 
2m dxj 



DjU, x'Ju = XjU, {xj^j 



■{xjDj + DjXj). 



Let X be a linear symplectic transformation xiv^v) = (^jO- The Segal formula (see 
e.g. the theorem 18.5.9 in [10]) asserts that there exists a unitary transformation M 
of L^(R'^), uniquely determined apart from a constant factor of modulus one, which 
is also an automorphism of ^(R") and ^'(R") such that, for all a G ^'(R^'^), 

(3.2) (a o = M*a'"M, 

providing the following commutative diagrams 



M 



M* 



and if a'" G C{L^{W)) 



M 
'2 



(aox)" 



M* 
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3.2. The metaplectic group and the generating functions. For a given x, how 
can we determine M ? We shall not need here the rich algebraic structure of the 
two-fold covering Mp{n) (the metaplectic group in which live the transformations 
M) of the symplectic group Sp{n). The following lemma is classical (and also easy 
to prove directly using the factorization of the lemma 12. ip and provides a simple 
expression for M when the transformation x has a generating function. 

Lemma 3.1. Let x = '^a,b,c be the symplectic mapping given by (12.41) . Then the 
Segal formula (13.21) holds with 

(3.3) {Mv){x)= /"e2^'^^("''')?7(r/)rfr/|det5|i/^ 



where S is given by (12.51) . 

3.3. Explicit expression for M. 

Lemma 3.2. Let x be the symplectic transformation o/M^ given by (I2.18p . Then 
the Segal formula (13.21) holds with M given by 

(3.4) {Mv){xuX2) = (AiA2)-'/'eM(^i^2)-i-c)xix2 

X [[ e2*"^"'''^''2t;(r/i,r/2)e2^'^(^^'^i''i+^2"'^2r?2)^^^^^2, 



(3.5) {Mv){xi, X2) = (AiA2)"'/'e''"<(^^^2)"'-^)"i"^(e'*"''"'^i^^tO(Ar'xi, A2 'xs). 

Proof. We apply the lemmas 13.11 and 12.41 along with the fact that the mapping 
Mp{n) 9 M X G Sp{n) is an homomorphism or more elementarily that (13. 2p 
implies for Xj ^ Sp{n), 

{a 0x2° XiT = M*{a o X2TM1 = M*M;a'"M2Mi. 

The factorization of the lemma 12.41 implies that 

(M^;)(x) = e^"<^^'^> / e2^'^<^^'''>e*^<^'''^>t?(r/)rfr/, 

which gives readily the formulas above. □ 

Summing-up, we have proven the following result. 

Theorem 3.3. Let q be the quadratic form on given by (I1.12p . We define the 
symplectic mapping x by (12.180 and the metaplectic mapping M by (13. 5p . We have 

(3.6) (g o x)iy, v) = lAvl + lAvl + + vl^ {^he fi] are given by i^^), 

(3.7) (g o xT = M*q'"M. 

We can also explicitly quantize the formulas of the lemma 12.61 to obtaiiil 



2 2' 
(3.8) = (^(Aicd - dX2^)x2 + AiD,,) + /i^ (^X2d~^D.,, + cAsXi' 

-|- (^{X2cd — dXi^)xi + A2-Dx-2 j + (^'^I'^^^-Dxi + cAiX2 



^Note that for a linear form L on M^", L""!,"' = (L^)™ 
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4. The Fock-Bargmann space and the anisotropic LLL 
4.1. Nonnegative quantization and entire functions. 

Definition 4.1. For X,Y E M^" we set 

(4.1) n(x,r) = e-fi^-^i'e-*-[^'^i, 

where [X,Y] = {aX,Y) is the symplectic form (12. ip . For v G L^(]R"), we define 

(4.2) {Wv){y,T]) = {v,ipy,,)LHR^), with ^y,,{x) = 2«/4e-(-2^)'e2-(--i)^. 
We define also 

(4.3) Ao ={u e /^^(ffiy"^) such that u = /(z)e"5l^l', z = r] + iy , f entire.} 

Proposition 4.2. The operator Uq with kernel Il{X,Y)is the orthogonal projec- 
tion in L^(M^") on Aq, which is a proper closed suhspace of L^(M^"), canonically 
isomorphic to L^(M"). We have 

(4.4) Ao = laJiW = L2(M2") n ker(a + |z), 

(4.5) W*W = ldi2(^,i-^ (reconstruction formula u{x) = / {Wu)(Y)ipY{x)dY ), 

(4.6) WW* = Ho, (W is an isomorphism from L^(R") onto Aq). 

Proof. These statements are classical (see e.g. [12]) ; however, since we shall need 
some extension of that proposition, it is useful to examine the proof. We note that 
e~^'^y^(Wv){y,T]) is the partial Fourier transform w.r.t. x of 

M" X M" 3 {x,y) ^ v{x)Tl^e-^^''-y^\ 

whose L^(M^")-norm is ||w||i2(]Rn) so that W is isometric from L^(M") into L^(M^"), 
thus with a closed range. As a result, we have W*W = IdL2(Kn), WW* is selfadjoint 
and such that WW*WW* = WW*: WW* is indeed the orthogonal projection on 
raniy (ranl^VT* C ran and Wu = WW*Wu). The straightforward computation 
of the kernel of WW* is left to the reader. Let us prove that Aq = ran W is indeed 
defined by 04.31) . For v G L'^iW), we have 

(4.7) {Wv){y,7])= I i;(x)2"/V"("-^)'e-'*"("-2)^da; 

= I t;(x)2"/^e-"(^-^+*'')'rfxe-^(3''+^')e-5(^+*^)' 
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and we see that Wv G L^{R'^"')nker{d+^z). Conversely, if $ G L^{R'^'^)nkeT{d+^z), 
we have ^) = e~5(^^+^^)/(^ + ix) with $ e L^(M^") and / entire. This gives 

= e-f «^+-^) j j e-f (,^+,^+2^,(^+«)-2,«+^.)) ^^^^ ^^^y^^ 

= e-f «'+^') e— fe'+'?')e"(''-*?')«+"-)/(^ + iy)dydri 

= e-^l^l' y I e~-\^\\<^fiC)dydv {( = V + ^y, z = ^ + tx) 



since / is entire. This imphes 
proposition is complete. 

Proposition 4.3. Defining 



= $ and $ G ran TV. The proof of the 

□ 



(4.8) = ker(9 + -z) n ^'(R^"), 

the operator W given by (14.21) can be extended as a continuous mapping from ^'(M") 
onto (the L^(M") dot-product is replaced by a bracket of (anti) duality). The 
operator U defined by its kernel U given by (14. ip defines a continuous mapping from 
=y(]R2") into itself and can be extended as a continuous mapping from =5^'(R^") onto 
J^. It verifies 

(4.9) = n, n,^ = id^ . 

Proof. As above we use that e~'^'^y'^{Wv){y,ri) is the partial Fourier transform w.r.t. 
X of the tempered distribution on R^^y 

^;(a;)2"/^e-"(^-^)^ 

Since e^'^'"^^ are in the space 0'm{R^^) of multipliers of ^(M^'^), that transformation 
is continuous and injective from ^'(M") into S^'{E?^). Replacing in (14. 7p the inte- 
grals by brackets of duality, we see that W{S^'{W^)) C . Conversely, if $ G J(f , 
the same calculations as above give (14.91) and (14. 8p . □ 
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For a Hamiltonian a defined on M^", for instance a bounded function on M^", we 
define a^i^"^ = W*aW : 

(multiplication by a) 



W 



W* 



L2(M«) . 

„Wick 

we note tliat a(x,^) > =^ a^"^*^ = W*aW > 0, as an operator. Tliere are many 
useful applications of the Wick quantization due to that non-negativity property, but 
for our purpose here, it will be more important to relate explicitely that quantization 
to the usual Weyl quantization (as given by (13. ip ) for quadratic forms. 

Lemma 4.4. Let q{X) = {QX,X) be a quadratic form on M^" (Q is a 2n x 2n 
symmetric matrix). Then we have 

(4.10) g^"^'= = g"' + — trace Q. 

An 

Let L[y, ri) = T- y — t-rjhea real linear form on M^"; then, for all $ G Aq, we have 



(4.11) // L{y,r,f\^{y,r,)\''dydr^> V lli^(M^^)- 



Proof. A straightforward computation shows that 

(4.12) gW'"'^ = (g * , where r(X) = 2"e-2'^l^l' (X G R'"). 

By Taylor's formula, we have (g * r)(X) = q{X) + 2"e-2^l^l' {QY, Y)dY, we can 
use the formula j^2^/H'^e~'^'^^^ dt = ^ to get the first result. For $ G Aq, we have 
= Wu with u G L'^iW) and thus 



//r2\Wick , //r2\?/i \ , trace(ij^) 11 

= {{L ) M,M)L2(]gn) = ((L ) M,M)i2(K„) H ||m 

and since L^L^ = (L^)"' for a linear form, we get since L is real-valued. 



2 t« 2 kP + I^P 2 

II-^^IIl2(M2") = 11-^"' ""1^2(18") "I II^IIl2(M2"); 

which implies (14.111) . □ 

N.B. The inequality (|4.11|) looks like an uncertainty principle related to the localization 
in fQj- functions of Aq. Moreover the equality (|4.1U|) provides a simple way to 
saturate approximately the inequality (|4.11|) : for instance if L{y,7]) = yi, we consider the 
sequence $e = Wue with Ue{x) = Lp{xi/ e)e~^/'^^{x'), ||v'IIl2(r) = ||'0||l2(R"-i) = 1, and 
we get, provided XLp{x) G L^'' 



j I yl\^e{y,r])\^dydii = j^xl\^[xi/e)\\-Hxi + i- = 0{e^) + ^. 
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4.2. The anisotropic LLL. Going back to the Gross-Pitaevskii energy (11. lip , with 
q given by fll.l3p . we see, using the theorem 13.31 and fl3.8p that, with u = Mv, 



2Egp{u) = {q'"u,u)L2(R2) + g j 

= {M*q^Mv,v)L2(^^2)+g j |(Mt;)(x)|^rfx 

= {iDl+fxiyl + Dl+fxlyl)v,v)L^(^R2)+g J \{Mv){x)\'dx 

= ( ((AiCfi — dX2^)x2 + AiD^-J^M + fil (^X2d~^Dx2 + cA2Xi)^m, m) 

+ {{{X2cd - dX^^)xi + X2Dx^)'^u, u) + (^^\{Xid^^ D^^ + cAiX2)^u, u) 

The question at hand is the determination of inf||tj|| ^=1 Egp{u), which is equal to 
miim^2=i Egp{Mv). Since /i^ = at e = (see (El])) and /i2 G [1,4] (see (O)), 
it is natural to modify our minimization problem, and in the {y, rj) coordinates, to 
restrict our attention to the Lowest Landau Level, i.e. the groundspace of Dy^+fj,ly2, 
that is the subspace of L^(M^) 

(4.13) LLLy = {vM ® 2i/V2^'e-'^'^^^'},,ei2(K) = ker(D,, - z/i22/2) n L\R^). 

If we want to stay in the physical coordinates (x, ^ ) we reach the following definition, 
obtained by using Segal's formula (13. 2p with M, x given in the lemma [3TT] so that 

LLL^ = M{LLLy). 

Proposition 4.5. Let q he the quadratic form on M"^ given by (I1.13p . We define the 
LLL as 

(4.14) LLL = (keiC) n L^(M2), with 

(4.15) £ = {X2cd - dXi^)xi + A2-DX2 - ifJ'2Xid^^Dx^ - 1^20X1X2 = V2 ~ ^1^21/2- 
The LLL is the subspace o/L^(]R^) of functions of type 

(2 2 \ 2 

- ^) + {P2X2)HI + ^)] j exp {-i^x^X2), 

where F is entire on C, and the parameters 7,/32,i^, a are given in the section^TT^ 
The real part of the phase of the Gaussian function multiplying -F(xi + 1(32X2) is a 
negative definite quadratic form when {uj^i^) 7^ (0,0). 

Proof. We have 

t^2y2 m 



iC = [^.2^id ^D^^ + /i2cAiX2 j +i (^A2-Dx2 - {d\ ^ - X2cd)x^ 

r,~ ( fJ'2^id~^di + 1X282 + 2z7r/i2cAiX2 + 27r{dX^^ — X2cd)xi 
2nx V / 

— {\^^2Xld~^^l + 1^X282 + i7r/i2cAiX2 + ii{dX^^ - X2cd)xx\ . 
TT \2 2 / 



ITT 

We set 



(4.17) ti = /i2^AiMxi, t2 = A2^X2, 
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and we get for z = ti + it2 
d 



dz 



d z — z z ~\~ z 

— + i7lH2C\i\2—^. h 7T{d\^^ - \2Cd)ljL2^id~^ ^ 

- + .vr- + z^-{l - 2\,\2C) = ^ + -vr- - z^-u a 



dz 



As a consequence, the LLL is the subspace of L^(C) of functions 
/(z)e"''^^^e''^^^'''' , with / holomorphic. 



We note that the real part of the exponent is 



^'-{ti + e2--{ti-ti)) 



71^2 



n ,2a; — z/^, n ,2a + z/^. 



2a 



2a 



and that 

2a - > [u, v) ^ (0, 0). 
Leaving the t-coordinates for the original x-coordinates, we get with / entire 

vr^2 



/(/i2 ^\ ^dxi + ^X2) exp 



2 



2a 



2a 



I.e. 



/(/.fg ^A-|^ ^dxi+iAg ^^2) exp 
and since 



vrAi2 



2a\l^2 



2..2^ + ^2V 2aA2 



exp (-«;r~r~^ a;iX2), 

2aAiA2/i2 



/i2Ai(i ^A2 ^ 



2-V2rf'ArV^' = 2-V27'4-'A^/i^" = 2-V27'4-'2/i27"'/52"V2 



^2Ai27-Ur'A2 ^A2~' = yU227"iA2"' = ^227"S/52(2/i2)~' = /32, 



-1.. .2/1-1x2,-2 



1.. .,2/|-l<- 



-1 ..-2 



7 

4/32 



2~V2A2~^ = 2"V2 

Tlll2^'^d 



1 7/^2 7/^2 



2^2 4 



2aAiA2/i2 2aAiA2 2a2 
we obtain 

f(^jj2^Xi^d[xi + i H2Xid~^ \2^ X2]) 

,2a 



7r/i2 

X exp 

that is, with F entire on 
(4.18) + i/32X2) exp 



2 j2/ 



Xid 



, 2a + z/^ 
2aAi ' 



exp [—1 



7lfi2l'^d 

2a\i\2fi2 



X1X2), 



'W2 



- ^) + (/^2X2)'(1 + — )J 1 exp (-il^xirra). 



.7rZ/2-y 



The proof of the proposition is complete. 



□ 
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Remark 4.6. We note that in the isotropic case z/ = 0, we have = 1,7 = 
4, recovering fll.lSp {f{xi + ix2)e~'^^^^~^^'^^) for a; = 1. On the other hand, the 
reader may have noticed that it seems difficult to guess the above definition without 
going through the exphcit computations on the diagonahzation of q of the previous 
sections. 

4.3. The energy in the anisotropic LLL. 

Lemma 4.7. The LLL is defined by the proposition "JTS and the Gross-Pitaevskii 
energy by (11. lip . For u G LLL, we have 



(4,19) EaA'") = i I ( , + 2l° + .^ ^'"' + a.-fJ+£' ^0 "^'1'''""'"^ 

g I 1 / \ |4 7 Ml Inn 1 \ 

+ -/ \u[xi,X2)\ dxidx2 + ^ P1P2 + -T-5- • 

Proof. In the LLL, one can simphfy the energy. We define 

A2 = M{r]2 - ifi2y2TM* = fi2 {Xid-^D,, + 0X^X2) + i {X2D,, - {dX^' - X2cd)xi) , 

Ai = M{r]i - i^iiViYM* = fii (Asc/'^D^.^ + 0X2X1) + i {{Xicd - c/A^^)x2 + AiD^J , 

which satisfy the canonical commutation relations: [/lj,y4*] = fij/n, while all other 
commutators vanish. We have proven that 

= AlA^ + A;A2 + ^^Ip^ = {ReA.f + (Im^O' + (Re + (Im 

ZTT 

and the LLL is defined by the equation A2U = 0. On the other hand, we have 

d^i^ KeAi — lmy42 = dX^^xi, dii2~^ Rey42 — ImAi = dA2^^a;2, 

and thus for u G LLL, since A2U = , using the commutation relations of the A/s, 
one gets 

d^X^^xl = d^fi^\ReAif + {{A2 - A;)/2ty + 2dfi^\Re Ai){A2 - A;)/2i 
= ciVr'(ReA)^ + £, 

and similarly, 

d^X^^xl = d^fi2\{A2 + A;)/2f + {ImA^f 

d^ 



(ImAi 



47r/i2 

As a result, we get on the LLL, 

fijXfxl + d'Xfxl = (ReAi)2 + (Im + + 



47r/i2 Aird"^ ' 



and = /i^A^^x? + d'Xfxl - ^ - fi^ + ^, so that 

2Egp{u) = 7: ( fiif3ixl + ^xl] \u{xi,X2)\'^dxidx2 
2 Jr2 V Pi J 



+g / \u{xi, X2)\^dxidx2 



2!r 4ir V Afi 
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for any u G LLL, that is, satisfying (14.161) . We note that 

'^^^ = — , I 2 ^^' (coefficient of xj) ^ = '^^'^'^'^ + ^ ^ (coefficient of xj) . 

Definition 4.8. For u G LLL (see the proposition^^, we define 

(4.20) £lll{u) = \ I {e^x\ + K\xl)\u{xi,X2)\^dxidx2 + ^ [ \u{xi,X2)\'^dxidx2, 

with 

,2 I ,,2VO,,2 I r-2\ ^, I o, ,2 I ,,2 



2 (a + 2u;2^z/2)(2z/2 + 52) « + 2cu2 + z/' . . ^ , 

a — + 2uj^ 2a 
We note that, from fl4.19p . 

(4.22) Egp{u) = -^^——Slll{u) + ^-^ ( f3iP2 ' ^ 



Remark 4.9. Since = z/^ + 4u;^, we see that 

X . 9 9x/ 2z/2 , „ (a + 2u;2 + z/2)(2j/2 + 
4.23 2Z/2 + 52 1 + = = ^ ^ > 2//^ 

and = z/ = 0. 

Remark 4.10. We stay away from the case where u = and shall always assume 
uj > 0. In the case u = 0, the quadratic part of the energy is diagonal and the LLL 
is, 

and we get a ID problem on the function vi. 

4.4. The (final) reduction to a simpler lowest Landau level. Given the fact 
that in fl4.16p . we can write F{xi + 1^2X2) as a holomorphic function times e~^^ , 
with 6 = 77rz/^/(8/32a), and that the energy Slll depends only on the modulus of u 
and not on its phase, it is equivalent to minimize Slll on the LLL or on the space 

f{xi + 1(32X2) exp {^-^ ^1 + {P2X2Y ^ , with / entire. 

A rescaling in xi and X2 yields the space of the introduction with 



1 I 1 1P2 

(4.24) u{xi,X2) = X -v{yi,y2), Vi = Xi. y2 = X2\ ' 



2 2p2 V 2 ' 

and, with Aq given by (14. 3p . the mapping LLL 9 m 1— *• u G Aq is bijective and 

isometric. With K,i,gi given in the definition 14.81 /52 in (12.121) . 7 in (I2.13p . we 
introduce 

(4.25) = IT^ 9o 



^2" 4/^2' 
and 

(4.26) E[v) = \l {eS\ + KS''^\v{y^.V2)\^ dy^dy2 + ^||^| 
Using the transformation (I4.24p . we have 

(4.27) Elll{u) = ^E{v), 

7 
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SO that, via the definition 14. 8^ we are indeed reduced to the minimization of fll.2ip 
in the space Aq (given in fll.221) ) under the constraint ||M||i2(]g2) = 1. We note also 
that the quantities 

2a 2p2 

a + 2uj'^ + i/^ ' 7 



(4.28) — — -, (factors of SLLLiu) in and E{v) in gSZ])), 



(4.29) and P2, — , + + (factors of k in (I12SD. of gi in g^S])), of g in , 

P2 2a 

are bounded and away from zero as long as u stays away from zero, a condition that 
we shall always assume, say < ujq < uj < 1. 

5. Weak anisotropy 

This section is devoted to the proof of Theorem II. 1[ We assume e < k, <^ e^^^. 
The isotropic case is recovered by assuming k = e. We first give some approximation 
results in subsection 15.11 and prove the theorem in subsection 15. 2[ 

We recall that the space Aq, the operator Ho, the energy E and the minimization 
problem I{e,K) are defined by ffT^ . (fT^ . ffT^ and fOij) . respectively. An 
important test function will be (11.281) . namely 

(5.1) u^(xi,X2) = e^(^'"l^l')e (v/r72;,r) , z = xi+ix2,. 

for T = tr + iTj = . 

5.1. Approximation results. 

Lemma 5.1. Let u{x) = /(xi + zx2)e~^'^''^ G L°°(M^), with f holomorphic. Assume 
< /3 < 1 and let p e C°'^(M2) be such that supp(p) G Bs the Euclidean ball of 
radius S > and of center 0. Define 

Then, for any r > 1, there exists a constant Cs,r > depending only on S and r 
such that, setting R = mm{Ri, R2) , we have. 



(5.3) IIIIo (pn) - pu\\j;^rm2) < Cs,r\\u\\L 



{RlR2y" 



RP ■ 

Proof. We first prove the lemma in the case (3 = Q. For this purpose, we write 

\Mpn)\ < I e-^^\^-y\'\u{y)\\pmdy- 



Young's inequality implies, for any r > 1 and any p, g > 1 such that 1/p + 1/g 
1 + 1/r, 



e 2i-''i 



II^^pIIl, < ||^^||l° 

LP 



LP 



||no(pn)||^, < 
Fixing q = r, hence p = 1, we find 

(5.4) ||no(/w)||^, < 2||M||i^ IIpII^, = 2||m||l^ {RiR2)-^^ 

This proves (15. 3p for /5 = 0. 
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Next, we assume P = 1. We use a Taylor expansion of p{y) = p{x + y — x) around 



x: 



p{y) = p{x) 



1 



VR1R2 Jo 



Ri i?i R2 R2 J \ Ri R2 



dt. 



We then notice that, ahhough n ^ Aq a priori, it belongs to J(f^ (see the proposition 
I4.3p and we have Ilo{u) = u since u G and u{x) = f{xi + 1X2) e'XJp{—7i\x\'^ /2) 
with / holomorphic. Hence, we have 



no(pu) - pu 



^S+1 



X 



yRlR; 



Vp 



2 JO 



xi , ^ yi - xi X2_ y2 - X2 \ ( yi - xi y2 - X2 



Ri Ri ' i?2 



R2 



- p{x) 



Ri R2 



dtdy, 



u 



^S+l ) 



where the set B^^^^ is 

(5.5) Sj;f = {(1/1,1/2) = (i?A,i?2t2), teBs+i}. 

We thus have, with R = min(i?i, R2), 

1 \y-x\ 



(5.6) |no(pn)-pn|<||Vp|h 



Hy)\ 



\/ R1R2 R 



dy 



+ |p(x)| 



|M(?/)|e-tl"-2^l'd2/. 



We bound the first term of the right-hand side of 05.61) using Young's inequality, 
while for the second term, we have, Vx G supp(p) C Bg^'^^, 



\u 



(B-|f^) 



{y)\e-'^\^-^\" dy < 



I II _2Lr2 



\x~yY 



dy 



II ,1 _ TT I~,2 II II C 

-fx 



where C is a universal constant. Hence, we have 



l|no(pM) -pwlLr < -^IIVpII^oo \y\e 



C 



\\U\Uoo—= 



R-Ri,-R2|l/r 



C 11 

+ ^l|M||L-|b|U-(i?li?2)^~^|5s| 



l/r 



This gives fl5.3p for /9 = 1. We then conclude by a real interpolation argument 
between C° and □ 
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A comment is in order here: we have chosen to state Lemma 15.11 with a general 
function p. However, since our aim is to apply the above result with the special case 



p{x) 



[I - \x 



2\l/2 



it is also possible to use explicitly this value of p in order to 



give a simpler proof of the above result. The method would then be to prove the 
estimate for r = +00 first, then for r = 1, and then use an interpolation argument 
between and For instance, the proof of the r = +00 case would go as follows: 



\Uo{fm){x) - p{x)u{x) 



Q-^\x-y\^+iTT(x2yi-y2Xi) 



{p{y)u{y) - p{x)u{y)) dy 



< 



< 



\u\\l^ 
\\u\\l° 



fl^-^l' \p{y)-p{x)\dy 



\x - y\ 



R 



R 
y\dy. 



dy 



The proof of the case r = 1 is slightly more involved, but is based on the same idea. 
We now prove 



Lemma 5.2. With the same hypotheses as in Lemma \5. il we have, for any s > 1, 



(5.7) 

and 

(5.8) 



1/2 



X^' \Ilo{pu) - pu\ ] < C5,s||m|U°o(M2)||p||cO,/3(m2) 



1 + RjS' 

rF~ 



1/2 



X2' |no(pM) - pu\ ] < Cs,s\\u\\l^{R2)\\p\\cO:P 



RP 



where Cs,s depends only on S and s. 

Proof. Here again, we first deal with the case P = 0. For this purpose, we write: 



(5.9) \x,\'\Uo{pu)\<r-' [ 



ki-2/ire-5'"-^l' 



u{y)\p{y)dy 



+ 2 



s-l 



\u{y)\p{y)dy, 



where we have used the inequality (a+b)'^ < 2*"^(a''+6*), valid for any a,b > 0, s > 1. 
The first line of (15. 9p is dealt with exactly as in the proof of Lemma 15.11 leading to 
(15 ■4p with r = 2, which reads here 



(5.10) 



|xi - yil'e 2!^ y\ \u{y)\p{y)dy 





1 IS --li:|2 


< \\u\\loo 


\x\ e 21-^ 


L2 





\P\\l^ 



< Cs\\u\\l^\\p\\l2, 



where Cg depends only on s. The second line of (15. 9 p is treated in the same way, 
but p{y) is replaced by \yi\'^p{y), that is, p{y) is replaced by Rl\yi\^p{y). Hence, we 
have 

(5.11) [ \y^\'e-^^\^~y\'\u{yMy)dy 



L2 
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Collecting 1^, CTID and (imi) . we find 

|||xi|^no(pn)||^, <C,(l + i?J^^)||n|Uoo||p||co|55r/l 

This proves (15. 7p for /9 = 0. 

Next, we consider the case (3 = 1. Here again, we use a Taylor expansion to obtain 
(15. 6p . This implies 



|no(pM) - pu\ < 2' 



-illVplh 



+2^ 



R 

-illVpll,, 



^S + 1 



A/-R1-R2 



1 



A/-R1-R5 



+ |xinp(x)| 



S+1 J 



|n(y)|e- 



where Bg]^^^ is defined by (15. 5p . We use Young's inequality again, finding 



\Ilo{fm) - pu\\\j^2 < 2^ 



-illVplh 



R 



-illVpIL 



R 



\y\e 21 



\ R1R2 

2s 



1/2 



1/2 



11/1 1 



JrC 



where C is a universal constant. Hence, 



\x,nUo{pu)-pu\\\^2 < CsJ^^il + RlS') 

ri 



This gives (15. 7p in the case (3 = 1. Here again, we conclude with a real interpolation 
argument. The proof of (15. 8p follows the same lines. □ 

5.2. Energy bounds. 

Proposition 5.3. Lei r G C \ M, let p e C^''^/^{R'^) be such that supp(p) C K for 
some compact set K , and J |pp = 1. Consider as defined by (ll.28p . and define 

(5.12) V = ||Ho(pn^)||^2(K2) no(pu^), 

where p is given by 



(5.13) p(x) 



1 ( X]_ X-A r> ^ ( ^90^\ 



Ri 



Then we have, with E{u) defined by (ll.2ip 

1 



(5.14) E{u) 



2geK 



x\^\p{x)\^ + ) + O I 



3\ 1 



for {e,Ke ^/'^) (0,0), where 7(r) is given by (ll.3ip . 

N.B. The L°° function pu^ does not belong to Aq since it is compactly supported 
and not identically 0; as a result, ||Ho(pUt)||l2 ^ and v makes sense. 
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Proof. First note that R = mm{Ri, R2) = R2, and that Lemma l5.ll with r = 2 
imphes 

(5.15) \\\Mpu^)\\l2 - Wpu^L'A < CR~''^ = C [^\ . 

We then apply Lemma [5.21 for s = l,f3 = 1/2, finding 



x'i \Uo{pUr)\^ - I xl\p\'^\u. 



1 -\- R 

< C{\\XiIloipUr)\\L2 + \\XlpUr\\L^)^Y/r 



We also compute 



Hence, we get 



(5.16) 



l + Rj 



.3\ 1/8 



A similar argument allows to show that 



K 



(5.17) - 



xl\Uo{pUr)\^ - I Xl\p\'^\u '2 



Turning to the last term of the energy, we apply Lemma 15.11 again, with r = 4, (3 
1/2, finding 



\Uo{pUr)f - / \pUr 



< 2 (||no(pn^)||j4 + ||p^^r|li4) llno(pn^) - pUrW^i 

< C\\pUr\\%iRlR2r'^'R-'^'. 



In addition, we have 



[ \pUr\' < \\Ur\\U [ \p\^=\\u4U{R,R2r' [ 



Hence, we obtain 



(5.1^ 



|no(pn^)|^ 



\ ^ / 



3\ 1/8 



Combining fl5.16p . 05.171) and fl5.18|] . we have 



E (UoipUr)) = EipUr) 

Hence, with the help of f l5.15p . we get 



1 + 



3\ 1/8' 
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Finally, we estimate the terms of i?(/7u^/||p'u^||i2): using real interpolation between 
C° and C^'^, we obtain 



(5.19) \\pu4l. 



\p{x)\'^\Ur{RlXi, R2X2)\'^dx 



/ki 



o 



\urr + o\ ( - 



3\ 1/8' 



Moreover, we have 



(5.20)- / x?|p|>,'2 



(5.21)y / 



(5.22) ^ / ipiv^r 



-Ro 



9 



2R1R2 



\Ur\^ + 

\uX + 



3\ 1/8' 



3\ 1/8' 



.3\ 1/8' 



xl\p{x)\'^dx, 
xl\p{x)\'^dx, 

\P\'- 



Thus, collecting fICTll . (^M>, flCTll and (Km . 



E{u) 



'■ r f 

-Rl xl\p{x)\^dx+ —Rl xl\p{x)\'^dx 



(^|^^|2)2 2i?ii?2 



25(oeK 



vr 



1+0 



25(o£:k 



7r 



3\ 1/8' 



.3\ 1/8' 



+0 I v^(^^ 

□ 

Proof of Theorem \ 1 . 1[ We first prove the lower bound in (11.331) : this is done by 
noticing that 

Ji^, 't) < 

where 



\u\' = iy 



J(5,fi:) = inf {^(m), uE L'^ {R\{l + \x\'^)dx)nL\R'^), I 

In addition, the minimizer of J(£:, k) may be explicitly computed (up to the multi- 
plication by a complex function of modulus one): 

„2 \ 1/2 
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with Ri,R2 defined by fl5.13p . Inserting (15.231) in the energy, one finds the lower 
bound of fll.33p . In addition, the inverted parabola fl5.23p is compactly supported, 
so it cannot be in Aq. Hence, the inequality is strict. 

In order to prove the upper bound, we apply Proposition 15.31 with 

1/2 



p{x) 



1 



\x\ 



TT 



and r = j. This corresponds to minimizing the leading order term of (15.141) with 
respect to r and p, with the constraint J |pp = 1. □ 

6. Strong anisotropy 

We give in this Section the proof of Theorem II. 2[ We deal here with the strongly 
asymmetric case that is, (ll.35p . which we recall here: 

(6.1) K>ei/3 

We first prove an upper bound for the energy in Subsection 16.11 then a lower bound 
in Subsection 16.21 and conlude the proof in Subsection 16.31 

6.1. Upper bound for the energy. 

Lemma 6.1. Assume that p G L^(M). Then the function 



(6.2) u{xi,X2) 

satisfies m G Aq. 
Proof. We first write 

u{xi,X2)e" 



21/4 



e 



:((.,-,,)^-2..x..)p(2/i)d2/i, 



2V4 



which is a holomorphic function of Xi + 1x2- In addition, we have 

\u{xi,X2) \ < ^e-5"2 
Hence, using Young's inequality, we get 



p * e' 



2 



[Xu 



1 



2l/4| 



P\\l2( 



IfIIl2(r2) < ^\\p\\m 

hence m G L^(M^). 

Lemma 6.2. Let p G C^(]R) have compact support with supp(p) C 
consider the function 



□ 

-T,T), and 



(6.3) 



Pit) 



p 



t 



Then, for any r > 1, there exists a constant Cr depending only on r such that the 
function u defined by (16.20 satisfies, for R> 1, 



(6.4) 



u 



'xi,X2) - 2^/V(a^i)e""''^+'"''''^' - z2^/^'^onV'r, V— ^i+^-^1^2 



X2P [xi)e 



< CrT^' 



^5/2-l/r- 
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Proof. We use a Taylor expansion of p (^) around that is, 
(6.5) p(|)=p(£l)+ip'(|)(j„-.0 



In addition we have 



2V4 
and 



Ry^ Vi?, 
Setting 

(6.6) v{xuX2) = u{xi,X2) - 2i/V(a;i)e""'^'^'"^'^'' - z2^/^X2p'(xi)e-"^2+i-^i^2^ 
we infer 




|y(xi,X2)| < 21/45/2 g"^"' / / yie-^^'(l-t) 



dtdyi 



Hence, using Jensen's inequahty, we see that there is a constant Cr depending only 
on r such that 



whence 



^ CrK^ [ [ t e-^^^^hfe-^y^l-t) [ l^.TR,TR){x,+ty,)dx,dtdx2dy, 



1 



= cJ^^^{2TR) [ I I e-^^^-hle-^^y\l-t)dtdx2dy, 

_ x-y/ Hp IIl°° fyn d 

which implies fl6.4p . □ 

Lemma 6.3. Under the same assumptions as Lemma lUTR let u he defined by (\6.2\\ . 
Then, there exists a constant Ct > depending only on T such that u satisfies 



(6.7) / xl u{xi, X2) - 2^/^p{xi)e-'''''^''''''''' - i2^/^X2p'(xi)e-^^'+'^^i^2 



dx 
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and 

2 



(6.8) / xl u(a;i,a;2) -2i/V(a;i)e-"^2+™2 _^2i/'^a;2p'(a;i)e-"^2+*-^'i 



dx 



„ lb"lli°°(R) 

Proof. Here again, we use the Taylor expansion f l6.5p . Hence, v being defined by 
(16. 6p . we have 

\xi\\v{xi, X2)\ < ^lyi^kile-^'^'y J y^e-^y^l - t)l^_TR,TR){xi + tyi)dtdyi 

- 2^^""'^"' / / ^1^^^'''^^ " + tyi\M-TR,TK){xi + tyi)dtdyi 



1 1 1 1 

\\P \\L°° ^ILr^ 
-I " " e 2 2 



/ / \yi\^e ^yh{l-t)l(^^TR,TR){xi+tyi)dtdyi. 
Jw Jo 

Hence, using Jensen's inequahty and arguing as in the proof of Lemma [6. 2 ^ we have 
where C is a universal constant. This implies (16.71) . A similar computation gives 



which proves (16.81) . □ 

6.2. Lower bound for the energy. We first recall an important result by Carlen 
[7] about wave functions in Ag (defined by (11.221) ): 

Lemma 6.4 (E. A. Carlen, [7j). For any u G Aq, Vn G L?' , and we have 
(6.9) / |V|n|p = 7r / 



Remark 6.5. The result of Carlen is actually much more general than the one we 
cite here, but the special case (16. 9p is the only thing we need. 

Lemma F6.4I implies the following decomposition of the energy in Aq: 

Lemma 6.6. Let u & Aq be such that \\u\\l2 = 1. Then, we have 



(6.10) / \dMf + '- I r.;H= + f / M'. 

Proof. We write 

(6.11) Eiu) = -^ + ^ + ^ [ xl\u\^ + '^l xl\u\^ + f [ \u\\ 

Sit 87r 2 J^2 2 J^2 2 J^2 

Hence, applying (16.90 . we find (16.100 . □ 
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Note that the first hne is easily seen to be bounded from below by the first 
eigenvalue of the corresponding harmonic oscillator, namely K^/(47r). Hence, f l6.10p 
readily implies 

(6.12) Eiu) > ^. 

This explains why we chose the constant |^ in the decomposition (16.111) : it is the 
constant which gives the highest lower bound in (I6.12p . 

6.3. Proof of Theorem II. 2L 

Step 1: upper bound for the energy. We pick a real- valued function p such that 
p e C2(M),supp(p) C (-T,T), [ p' = l, 



and define u by (16.21) . where p is defined by (16. 3p . with 

(6.13) R = e-'^l^. 

Hence, setting v = -n-jr — u, we know by Lemma 16.11 that w is a test function for 
/(e, k). Hence, 

(6.14) I{e,K)<E{v). 
Next, we set 

= 2^/V(a;i)e""'^'+^"'^^'^^ +z2i/^X2p'(xi)e-"^2+i-^i^2^ 
and point out that, applying Lemma [6.21 with r = 2, 



Mh = \\Vl\\l2 



+ O (5^/3^ = 1 + 2^/2 [ |p'(xi)|2 / x^e-2-^'dx2 + 0(5^/3^ 

Jr Jr 

= 1 + Ce'/' [ p'' + 0{e'/'), 
Jr 

where we have used that the two terms defining vi are orthogonal to each other. 
Hence, 

\\uh. = 1 + O {e'/') , 

where the term O {s^^^) depends only on ||p'||i2, and T. According to (I6.14p 

and the definition of f , we thus have 

(6.15) I{e,K)<E{u)[l + 0{e^/'')], 

where the term O {s'^^^) is independent of n. We now compute the energy of u: 
applying Lemma [6.31 we have 



2| |2 / 2| I 



< Ce^/^ (||xin||i2 + ||xit;i||L2) < Ce'^^^ (2||xifi||L2 + 
Moreover, we have, since p is real-valued, 

/ xl\vi\'^dx= I xlp{xiYdxi + ^ I xlp{xifdxi = e'^^^ ! t^p{tfdt + O {!) . 
Jr2 Jr 4:71 Jjj Jjg 



Hence, we have 



(6.16) / xl\u\^ = e-^/^ [ t^p{tfdt + O {1] 

Jr2 Jr 
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The same kind of argument allows us to prove that 



(6.17) 



2| |2 
X2\U\ 



Next, we apply Lemma F6.2I with r = 4: 



\u\ 



\Vl\ 



< 2||n - v4l^ {\\u\\l, + ||t;i||i4) < Ce^'^ {\\u\\l, + |K;i||i4) 



Moreover, we have < II'^iIIl-i + Ce^^^, hence 



\u\ 



Li- 



We also have 



.2/3 



Hence, we obtain 
(6.18) 



\u\ 



An 



647r2 



,'4 



£2/3 f p{t)^dt + 0{e''). 
Jr 



Collecting (l616ll . (ICTl) and (16381) . we thus have 

E{u) = f + {k^s'/') + e'/' ( I \epmt + I / p{tfdi\ + O (.2) 



Recalling (16.151) . this implies 



2 



£2/3 - 2 _ 

As a conclusion, we have 
lim sup 



r2/3 



9g 



< =- / t'p{tydt + ^ / p{tydt 



for any real-valued p G C2(M) having compact support, and such that ||p||l2 = 1. 
density argument allows to prove that 



lim sup 



I(e, k) - ^ 



£2/3 



where J is defined by (ll.37p . Thus, we get 



,2/3 



J + c[ e 



A/3 



K 



with lim c(t, s) = 0. 

(t.s)-.{0,0) ^ ' 
t,s>0 
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Step 2: convergence of minimizers. Let m be a minimizer of /(e, n). Then, according 
to the first step, we have 

E{u)<^ + Je'/' + e'/'c(e,— 

with hm c(t, s) = 0. Hence, applying Lemma [6.61 we obtain 

(t,s)-.{0,0) \ ' ' 7 r-jr J o 

t,s>0 



(6.19) 

+ } 

We set 
(6.20) 



2 V47r2 



|<92|^^| 



2 1 |2 



2 

in I l|2 , ^ / 2| |2 , 90 

\di\u\\ + — / x^\u\ + — 



U|4<^ + J^2/3^^2/3 / 



47r 



,1/3 



1 / Xi 

U I ^77:r,X2 



/3 



SO that ||f ||l2 = \\u\\l2 = 1, t> > 0, and (16.191) becomes 

^ ' |2 I / 2 2 

|c'2f| + / X2V 

2/3 



+ 



|<9lt' 



|2 , ^ 



This imphes that 
(6.22) 



\d2v\'^ / xy <c 



where C does not depend on (e, k). Moreover, since the first eigenvalue of the 
operator — + x\ is equal to l/(27r), (I6.2ip implies that 



(6.23) 



where C does not depend on (e, k). Hence, up to extracting a subsequence, v 
converges weakly in and weakly in to some limit vq > {]. Using fl6.22p and 
(16.231) . we see that 



x?v^ < C, 



hence v converges strongly in L^. Since in addition converges weakly in L^, we 
have: 



(6.24) 



V — > vq strongly in L 

(£,£1/3^-1)^(0,0) 



XiV 

V 



xiVq weakly in L^(]R^), 



(£,ei/3K-i)^(0,0) 

— > Vq weakly in L^(M^), 

(£,£1/3^-1)^(0,0) 

d2V — > d2Vo weakly in L^(M^ 

(£,£1/3k-1)-(0,0) 
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Hence, we may pass to the liminf in the two first terms of (16.211) . getting 
(6.25) 

t\ / \d2Vof+[ xy,< liminf ( -^f \d2vf + [ xlvA < ^. 

We use that the first eigenvalue of the operator —^-^ + on L^(M) is equal to 
l/(27r), is simple, with an eigenvector equal to 2^/^ exp(— 71X2). Thus, 

(6.26) ^;o(xl,X2)=e(a:l)2l/V'^"^ 
with ^ > 0. Next, flOT]) and flelMj) also imply 

(6.27) I I x\vl + f / vt< liminf (\ [ + f [ < J. 



Using (16.261) . we infer 

^ / ^'te + ^ I e<j. 



,2 



Hence, recalling that, in view of (16.241) and (16.261) . we have J = 1, the definition of 
J implies that ^ is the unique non- negative minimizer of (11.371) . This proves (11.381) . 
with strong convergence in and weak convergence in L^. Moreover, using (I6.27P 
again and the fact that ^ is a minimizer of (11.371) . we have 

lim ( [ xlivl - v"") +go [ U - v^)] = 0. 

(,,,1/3,-1)^(0,0) \u ° ^ ° V 

Next, using the explicit formula giving vq, a simple computation gives 

xliv' - v',) + goiv^ - v',) > g [ {v^ - v^.f 

hence v"^ converges to Vq strongly in L^(]R^). Thus, 

4 I i 
V — y Wo- 

The space L^(]R^) being uniformly convex, this implies strong convergence in L^, 
hence ( lL38ll . 

Step 3: lower bound for the energy. Using Lemma 16.61 we have 

E{u)>^ + '-^([ xy + go [ 
4vr 2 \J^2 Jr2 

In addition, we already proved (ll.38p . which implies 

l/_22,fi'o/ 4 1/ 2 2 . 9q 



^ JlR2 / J]R2 Z J^2 Z Jk2 

which implies the lower bound for the energy. □ 



7.1. Glossary. 



7. Appendix 



38 



A. AFTALION, X. BLANC, AND N. LERNER 



7.1.1. The harmonic oscillator. The operator 

(7.1) E + X]x]r = E ^(Dl + A?x|), A, > 0, D., = 
has a discrete spectrum 

(7.2) ^^ -^{E "A} 

and its ground state is one-dimensional generated by the Gaussian function 

(7.3) <^a(x) = 2'^/^ II AfV'^^^"?. 



(ai,...,a„)GN" 



7.1.2. Degenerate harmonic oscillator. Let r e {1, . . . ,n}. Using the identity 

(7.4) {Hru,u) = E ((Dl, + X'x^)u,u) = E IK^-, - ^A,x,)u||i. + ^Iklli^, 

we can define the ground state Ej. of the operator Hj. as 

(7.5) Er = L2(R") ni<,<, ker(D,^, - tXjXj) 



The bottom of the spectrum of nHr is | J2i<j<r -^i- 



7.2. Notations for the calculations of section 12.31 

(7.6) i/2 + a;2<l, z^^ + + = 1, 



(7.7) a = \/z^4 + 4^2 = ^4^2 + (i _ ^2 _ ^2)2 (if 1/ = o,a = 2u;). 

2 2 (1+Cc.2)2_^2 (i_^2)2_^4 2l/V + 6^ 
7.8 ;,2^ l+c^2_^^ ^ J J ^ 

(7.9) /i| = l + w^ + a (if i/ = 0,;U2 = 1 + w). 

Remark 7.1. If i/ = 0, /xi = 0(^2) and if / 0, /ii = 0(e). Moreover, for z^2 _)_ -|^^ 
G [1, 4] and for ly"^ + u}'^ = 1, ^2 g [2, 4]: we have indeed 

(7.10) 1 < l+w2 + (i/^ + 4a;^)^/2 <4 

since + 10cj2 < (1 _ ^^2^2 _^ iq^2 = 3^2 _^ ]^ _^ ^4 < 9 _^ ^4^ implying (3 - uj"^)"^ > u"^ + ^u? 
and (171^ . If z>2_^^2 ^ we have (l-u;^)^ = ^^"^ < i^^+4w2 ^ 2 < l+cj2 + (i.4_^4^2)i/2_ 



ANISOTROPIC BOSE-EINSTEIN CONDENSATES 



39 



We define the following set of parameters, 

r^^^ \ a ^^^1 a - 2u? - 2^2 4 2 4 2 2 2 

(7.11) Pi = 7, ^ = since (a - 2i^2)^ - j/^ = 4i^2 +4c^ - 4a;2Q = 4Lj2Atf 

a — 2uj^ + i/^ 2a;/ii 

(7.12) 132 = ^^(^^ _ a + 2u; since (a + 21^^)^ - = 40;^ + 4a;* +4a;2o = 40;^/^! 

a + 2uj^ + I/"' 2ci;//2 
2a 

7.13 ^ = _, 

a; 



(7.14) At 

(7.15) \l 



2 /ii 1 1 a — 2u? + 1^^ 

~ /ui + /3i/32/i2 ~ 1 + ^i^ ~ 1 + "+2'^^-'-^ ~ 2^ 

' 111 a—2ui^+u^ 

2 _ H2 _ 1 _ 1 _ a + 2u;2 + 



/X2 a+2aj^+!^ 

2 , ,,2n2 ^, ,4 



(7.16) and we have A? + Ai = l + -, A^A^ = 



!i2 , \2 



(7.17) d=^, c=^i^ so that cd=Ml^=^ 
^ ^ 2 ' 2A1A2 



We have also 



2/ii a - 2u;2 + a - 2^^ + z^^ 



2 



7/9i u;7 2a 

2/12 _ a + 2u;2 + i/^ _ a + + ^'^ _ -^2 
7/32 0^7 2a ^' 



and 



2Ai Va - 2u;2 + 



2 



(l + z.2a-i)2-i2V2ai/2 



- 2u;2 + u'^^/a + 2u? - v"^ 
\Ja + 2^2 _ j^2 



V^a2 - (2a;2 - z.2)5 



= (1 + z.2«-i)2-i/2c,i/2^ct + 2u;2 - z.2(2a;)-i(2z/2 + e2)-l/2_ 
Moreover, we have 



(7.18) CA2 = 2-3/2(«l/2 + ^,2^-1/2)^-1 



a + 2a;2 — 
2i/2 + £2 



A2d-i = ^ = 2-3/2(aV2 + i.2„-i/2)^-iy 



(ifi' = 0,cA2 


= 2-1/2, 


'a + 2u;2 - z/2 


2a; 


2zy2 + £2 


a + z/2 



A2d-1 = 2-V2(«l/2 + .^a"V2)^/^±|^^(« + ,2)-l^ 



(7.19) Asd^i = (2a)"i/y^^^J^^^ (if . = CA^rf-^ = 2-1/^(1 _ c.)-!/^), 

(7.20) cAi = ^i-t^ = (l + a-li.2)2-l/2„l/2(^ + 2w2^i.2)-l/2 (if ^ = 0, cAi = 2-V2(i + ^)-l/2)_ 

2A2 

Aid-i = Aic(cd)-i = (1 + a-V2)2-i/2aV2(a + 2a;2 + j,2)-i/22^(„ + ^2)-i^ 

(7.21) Aid-i = 2i/2a-^/2^(a + 2a;2 + ^2^-1/2 ^ ^ ^ 2^1/2(1 + ^)-i/2)^ 
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Xicd = 2-3/2(a + u^)uj-^a-^'\a - 2u? + v'fl'' (if = 0, Aicd = 2-1/2(1 - u^fl''), 
d 7A1 



au:-\a - 2u? + i.2)i/22-i/2«-i/2 = 2'^!'' a'l''u,-\a - 2u? + z.2)i/2^ 



X^cd -± = {a-2cO^ + Z.2)l/2 (2-3/2 ^ z.2)^-l^-l/2 _ 2-V2^1/2^-lN 
A2 

X,cd - ^ = 2-'/^u;-'a-'/\a - 2lo^ + z.2)i/2(c, + ^2 _ 2^)^ 
A2 



A,cd - f = -2-3/2a;-la"l/2(^ _ 2^2 ^ ^2)l/2(^ _ ^2)^ 

A2 



(7.22) 

(if y = Q, Xicd - ^ = -2-1/2(1 - t^)-i/2), 



(7.23) 



Ai 



-1/2 



a 



-l/2(c, _ 20,2 + ^2)1/2 (if ^ = 0, Ai = 2-1/2(1 _ ^)l/2)^ 



A2C(i - -f- = A^ iA2(AiC(i - -^) 
Ai A2 

= -2-3/2^-ia-i/2(a - 2a;2 + z.2)V2(„ _ ^2)(^ ^ 2a;2 ^ ^2)i/2(^ _ 2^2 ^ ^2^-1/2 

= -2-3/2^-i«-i/2(„ _ ^2)(^ ^ 2a;2 + z.2)l/2^ 

(7.24) X2cd-4- = -2-3/2^-^a-l/2(a-i/2)(Q, + 2tj2 + i.2)l/2 (jj ^ = q, Ajcd - — = -2-l/2(l + a;)l/2) 
Ai Ai 



(7.25) 
(7.26) 



-1/2 



-i/2(c, + 20,2 + ^2)1/2 (if ^ = 0, A2 = 2-1/2(1 + u;)i/2), 



2a 



a 



+ 2u;2 + z/2 



7/ii 4au;(2i/2 + ^2 



2/3i a - 1/2 + 2^2 



7.3. Some calculations. 

7.3.1. Proof of the lemma \KR. We have to calculate 

/ Ai 

n 1-1-7/2 n 

Q = x*Qx = X 





-Z/2 















1+1/2 


UJ 












1 





V 


—UJ 








1/ 





A2 

I J- — A2Cd CA2 

Xicd '0 cAi/ 



' 




-cAioj 



(l-Hl/2)A2-Hfi^-A2Cd<..; „Il±i^!l^+^cA2 



i.^A2 + T^-A2cd 



y -ijAi + y^-Aicd 

/ Ai 







LJA2 





-cA2 












An ^-A2Cd 



A 

A2 



CA2 




X 



— Xicd\ 



cAi / 

/(l-i/2)Ai-^+Aicda; 





y -LjAi + y^-Aicd 



'.^A2 + Y^-A2cd 









Ml^-eA,.\ 






d \ 
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We get easily qi2 = ^13 = = ^24 = 934- To prove that the symmetric matrix Q is 
diagonal, it is thus sufficient to prove that = = ^23- We have 

A^ 2\ \2 , Ai cdXi 2 2\2t 

gi4 = -{1 — u ) — uJcAi + cj-; — I ^ A^cuj — c A^a 



A? 
d 

A? 
d 



_J.(^l-iy^)-^cAi+UJ-^ + —^ 

d A2 A2 

2 1 cd'^ 2 79 



A2A1 A2A 



A? 




du^ 


A? 






A? 




du^ 


A? 







z/^ — a — z/^ + a 



(a;2 + u'^a) 



[a + 1/ J a 
2uj uo 

2\2-, 



(a + V 
4cu2 



2\2 



a + 1/ 



fz/^ + 4a;2 + z/2q,) 



fz/^ + 4a;2 _^ j^2^N 



2 + z/4 + 2az/2 



(z/^ + 4cj2 + + 2az/2) 



2 4 
(2z/4 + 8cj2 + 2z/2a) {2u^ + 4cu2 + 2az/2) 



4 



0, qed. 



Moreover we have 
A2 



^~,3 = _2l(l + ^2)^^^;^2 '^A2 , crfA2 , ,2,,„ x2„2. 



+ 



+ Aocuc — AoC d 



A, ' Ai ' -^^^^ ^^2^ 



A| 
ci 

A| 

d 



1 - z/^ + 2u;cd 
1 



ujd 



z/^ + a + z/^ — a + 



A2A1 A2A1 



2ci; uj 



(tt + Z/' 



2\2 



4aj2 



0, from the previous computation. 



We know now that Q is indeed diagonal. We calculate 



g44 

944 



A2(l-z/2) ^ 2cA?.. 2,2 _ A? 



rf2 



+ 



a;2d2 



^ c/2 

^z^^ + 4cj2 + z/^ + 2az/2 



1 -z/^+2cj c(i+c^(i^ 
A? 



Af 

rf2 



4cu2 



a;2(i2 



2u;2 + ct^2 



(z/4 + az/2) 



Since = 2d 2 = = 2/ io°2i 2m we have 



944 



a(a + 2tu2 + z^2) L 
Analogously, we have 

A2(l + z/2) 2cAlu 



Auj^ + 2auj^ + z/^ + az/^ 



933 
933 



A2 



rf2 



fc A2 - ^ 



l+u^-2ujcd+c'^d'^ 



^2 + 2au;2 + az/2 
a2 + 2aix'2 + az/2 



A2^ 



1. 



_2 

d2 



W2rf2 



2 9 z/^ + 4cj^ + z/* + 2az/ 



2^^ A2 



UJ^d^ 



1+z/ -a-z/N 
2u;^ - au;^ + 



a + V 



2\2- 



z/4 + az/2) 



Since ;3|2 = = ^ = a2(^_2^+.^) > have 



1 



933 



a{a — 2a;2 + z/2) . 



4u;^ - 2au;^ + z/^ + az^^ 



2au;2 + az/2 



a2 — 2aa;2 + q;z/2 



1. 
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We calculate 



udXi 

2ujd 
A1A2 



+ 2\lcduj + T2 - 2 

^2 



+ 2cduo + 



A1A2 A1A2 

2 



gii = A? (1 - - 2a + a + + ^ - 2 



a + a (a + i^'^Y 



2uJ uj 



+ 



A? 



n N 2 , 2 2 2 , (" + '^^)^ 
(1 — a)uj + a — a — au H 



a-2uj^ + 

2auj^ 
a-2uj^ + 



2aa;2 



oj"^ — au^ — + 



2uj'^ — auj^ av^ H 

2 2 



More calculations: 



= -8cu^ - 2cu2z/4 + a{2uj^ + 2cu2) 



which is equal to 

2acj2(l + tu^ - a) = a(2cj^ + 2cj2) - 202^^2 = a{2uj^ + 2cj2) - 2cj2(z/4 + 4^^), 

proving thus that gn = 1 + 0;^ — a. The previous calculations and (12. 8p give (/?g22 = 
1 + uj"^ + a, completing the proof of the lemma. 

7.3.2. On the symplectic relationships in Lemma \2.b\ The reader is invited to check 
the following formula^, with the notations of lemma 12.61 



{^1 



{6 



2. >2,6 + , 2, 



,a — Q, _|_ jy2 ^ 



2oj 



2uj 



2uj 



,a — V 



2oj 



2uj 
,q; — z/2 

^ 2oj 



=0, 



as well as 



. a - 2u;2 _^ ^2 ^ _^ 2ci;2 _ //2 ..^ 
( ) ( ^2 ^ ) «^ 



2"^£/i2^cu~^(a2 - (2 



2^2^1/2 



2a ^ 2a/i2 



and 



( 



a + 2cj2 + z/2 
2^^ 



)^/^2V2^( 



1 + u;^ + a 
a(a + 2u;2 + i/2] 



,1/2 



-^ = (l + o;2 + a)i/2 = /X2. 



This is indeed double-checking since those formulas are proven in section 2. 
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